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C\| i ABSTRACT 

We re-examine vacuum oscillation solutions of the solar neutrino problem taking (a) 
the results on total rates, electron energy spectrum, and the seasonal variations from 
the 708 day SuperKamiokande data and (b) those on total rates from the Chlorine and 
Gallium experiments. Best fit values for the mixing angle and mass splitting are found 
for oscillations to sequential and sterile neutrinos and the 90% CL. allowed regions are 
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1 Introduction 

Within the past year much has taken place in the arena of neutrino physics. The an- 
nouncement by SuperKamiokande (SK) of the strong evidence in support of a non-zero 
neutrino mass and oscillation in their atmospheric v data has catalyzed explorations 
of the subject from many different angles. The far-reaching impact of a neutrino mass 
on physics beyond the Standard Model and in astrophysics and cosmology is being 
vigorously examined. Additional details of the neutrino mass and the determination 
of the complete mass spectrum are therefore awaited with interest. 

The results on the arrival rates of solar neutrinos || from the Chlorine, Kamiokande 
and Gallium experiments were themselves strong indications of a non-zero neutrino 
mass. Recent high statistics results on solar neutrino rates, the scattered electron 
spectrum, and the seasonal variation of the rates from SuperKamiokande have added 
a new dimension to this effort @]. In this work we consider the vacuum oscillation 
scenario in the light of this body of data from several directions Q We utilize the latest 
708 day SK data |J and for the neutrino fluxes use the BP98 solar model which 
incorporates the INT normalisation. We consider the two alternatives of oscillation of 
the v e to a sequential (z/ M or i/ T ) or a sterile neutrino. 

2 Oscillation Probability 

In this work we restrict ourselves to the simplest case of mixing between two neutrino 
flavors. The u e survival probability used by us is 



where ■§ denotes the mixing angle in vacuum and A is the vacuum oscillation wavelength 
for neutrinos of energy E v given by AtcE u / Am 2 , in which Am 2 = \m\ — m\ | is the mass 
4 Prior to SK, the solar neutrino data have been examined in terms of vacuum oscillation by many 





authors Q. 
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square splitting. Here r is the distance of the point of neutrino production from the 
center of the sun and R(t) is the sun-earth distance given by, 



1 — e cos(27T — 

V rp 



(2) 



R(t) = Ro 

Here, Ro = 1.49 x 10 13 cm is the mean Sun-Earth distance and e = 0.0167 is the 
ellipticity of the earth's orbit, t is the time of the year at which the solar neutrino flux 
is measured and T is 1 year. 



3 Observed rates and neutrino oscillations 

After the declaration of the SK results, analysis of the total rates of all the experiments 
in terms of vacuum oscillation has been considered in || using the 504 day data. In 
this section we update this analysis by using the data collected over 708 days 0. The 
data that we use in our analysis for the total rates are given in Table 1. For the Ga 
experiments we take the weighted average of the SAGE and Gallex results. Because 
SK has better statistics we do not include the Kamiokande results. We fit the total 
rates from the three experiments using a standard x 2 -fitting procedure || . 
The definition of x 2 used by us is, 

X 2 = E (K th - F? p ) (^ 2 ) (Ff - ]<?*) (3) 

i,j=l,3 

Here F" = rB p 98 where a is th (for the theoretical prediction) or exp (for the experi- 

i 

mental value) and Tj is the total rate in the zth experiment. F^ xp is taken from Table 1. 
The error matrix contains the experimental errors, the theoretical errors and their 
correlations. For evaluating the error matrix we use the procedure described in [|Kj . In 
the presence of neutrino conversions, the detection rate on earth for the radiochemical 
experiments 37 CI and 71 Ga is predicted to be: 

(T)f = E / <j>k{E v )a{E v ) < P ee (K, r, R(t)) > dE u (4) 

h. J Eth 



where cr(E u ) is the cross-section for neutrino capture [11 and <fik(E u ) is the neutrino 



spectrum for the kth source []iT . The sum is over all the individual neutrino sources. 
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< P ee (E u ,r, R(t)) > is the average neutrino survival probability over one year for the 
time averaged total rates and further where a weighted sum over the production point 
in the sun has been carried out for each source. The theoretical prediction according 
to the BP98 solar model, T^ P98 , is obtained by setting the survival probability as 1 in 
the above. For SK, in the case of oscillation to sequential neutrinos one has to take 
into account the possible contributions from the or v T channel, 

(T)f = / dE A dE T p(E A ,E T ) dE v (j)i(E v ) 

JE A+U JO JE U - 



< P ee {E v ,r,R{t)) > ^+ < P efl (E v ,r,R(t)) > ^ 



dEr ' — ef .v-^- 5 --w/ - dEj , ( 5 ) 

The second term in the bracket is absent if oscillation to sterile neutrinos is under 
consideration. and E A denote the true and apparent electron energies respectively. 
p{Ea, Et) is the energy resolution function for which we use the expression given in 
[Q. E& th is 6.5 MeV for the calculation of the total rate at SK. is the differential 
cross-section for the production of an electron with true relativistic energy Et and can 
be calculated from standard electroweak theory. Below we summarize our results for 
the total rates. For the sequential neutrino case the best-fit values of Am 2 , sin 2 (2$) 
and xLin obtained are 

• Am 2 = 9.5 x KT 11 eV 2 , sin 2 (2tf) = 0.87, xLn = °- 52 

For 1 degree of freedom (3 experimental data points - 2 parameters) this solution is 
allowed at 47.08% CL. The best-fit point that we find is slightly different from that in 
[HI. We attribute this to the fact that in the above analysis we have not included the 



detector efficiencies. Incorporating these fl3|l , we get 

• Am 2 = 7.8 x 10- 11 eV 2 , sin 2 (2tf) = 0.74, xLn = 2 -54 

These values are in agreement with and with the 708 day data the quality of the fit 
is actually better. In fig. 1 we show the 90% CL. (x 2 < Xmm+ 4.61) contours for the 
vacuum oscillation solution. Apart from the global minimum there are several local 
minima as shown in the figure. The best fit for the sterile neutrino alternative gives 
X 2 = 5.74, which, for 1 degree of freedom, is ruled out at 98.3% CL. 
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4 Observed spectrum and neutrino oscillations 



In addition to the total rates, SK has provided the number of events in 17 electron 
recoil energy bins of width 0.5 MeV in the range 5.5 MeV to 14 MeV and an 18th bin 
which covers the events in the range 14 to 20 MeV 0. In this section we analyze the 
spectral data in the light of neutrino oscillations in vacuum. The definition of x 2 used 
by us is, 

X 2 = E (Xn^ - RT P ) ^ 2 (x n Rf - RT P ) (6) 

i,i=l,18 

where X n allows for an arbitrary normalisation of the 8 B flux with respect to the SSM 
prediction and Rf = Sf / S^ P9& with a being th or exp as before and Si standing for the 
number of neutrinos in the zth energy bin. The theoretical prediction is given by eq. 
(HP but the integration over the apparent (i.e., measured) energy will now be over each 
bin. Following SK we include the statistical error, the uncorrelated systematic errors 
and the energy-bin-correlated experimental errors as well as those from the calculation 
of the expected spectrum [ 14| . Thus the error matrix used by us is 



Cjj — fiij{Pi lS tat ^ijUncorr) a i,expO 'j,exp + Oi^al&j^cal (7) 

Our results are given below: 

(a) Oscillation to a sequential neutrino 

The best-fit values of parameters and Xmin are 

• Am 2 = 4.15 x 10^ 10 eV 2 , sin 2 (2^) = 0.89, X n = 0.75, xLn = 13.07. 

This solution is allowed at 59.68% OL. The allowed values of the parameters Am 2 and 
sin 2 (2$) at 90% OL. ( X 2 < xLn + 6 - 25 ) are shown in fi S- 2a - The best-fit values of 
Am 2 and sin 2 (2??) that we get are in agreement with those obtained in for vacuum 
oscillations. 

(b) Oscillation to a sterile neutrino 
For this case we obtain, 

• Am 2 = 4.16 x 10- 10 eV 2 , sin 2 (2$) = 0.77, X n = 0.76, x min = 12 - 9 



This corresponds to a goodness of fit of 61%. Thus, as far as the electron recoil 
spectrum data is considered the sterile neutrino alternative gives a marginally better 
fit. In fig. 2b we show the allowed values of parameters for this case at 90% C.L.f] 

5 Combined fits to rates and spectrum 

In this section we present the results of the combined fit of the total rates and the 
spectrum data treating the rates and the electron spectrum data as independent ||. 
Thus our definition of \ 2 is> 

x 2 = E - Fr) (i? - f;*v) 

+ E {XnRf - RT P ) °i? {X n Rf - (8) 

where the first term on the r.h.s is from eq. (f|) and the second from eq. (||). We allow 
the normalisation of the 8 B flux to vary as a free parameter. The Xmin an d the best-fit 
values we obtain are: 

(a) Sequential neutrino case 

• Am 2 = 9.47 x KT 11 eV 2 , sin 2 (2tf) = 0.80, X n = 0.71, xLin = 18 - 66 

For 18 (= 21 - 3) degrees of freedom this solution is allowed at 41.3% C.L.. In fig. 3a 
we show the 90% CL. allowed regions for the combined analysis of rate and spectrum. 

(b) Sterile neutrino case 

• Am 2 = 8.86 x lO" 11 eV 2 , sin 2 (2tf) = 0.88, X n = 0.82, xLn = 18.54 

5 In the analysis of the spectrum data we have not included the efficiencies as they are available 
in |i~3[ ] as functions of true energy for particular values of the threshold apparent energy. Here, we 
perform bin by bin integration over the apparent energy and there are different thresholds for each 
bin for which the appropriate efficiencies are unavailable. We do not have any reason to believe that 
the inclusion of the efficiencies will change the best-fit values dramatically - see e.g. the fits in Sec. 3 
- though the quality of the fit may be affected. 
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The goodness of fit in this case is 42.06%. In fig. 3b we show the 90% CL. allowed 
region for oscillations to a sterile neutrino from the combined analysis of rates and 
spectrum. We find that only two of the six regions allowed for the sequential case 
remain for the sterile alternative while for the others either the \ 2 is too high or they 
merge with these two. 



6 Seasonal variation and neutrino oscillations 

In [H the SK collaboration has also presented the preliminary seasonal data which 
shows a variation apart from that expected from the 1 / R 2 {t) dependence of the neutrino 



fluxes. If this is confirmed then it can help to distinguish between the MSW [L5|] and 
the vacuum oscillation alternatives. To analyze the seasonal data we define our \ 2 as 

2 



N th N exp 



X 2 = E — - 2 ' 7 (9) 



M=l,8 ° 



Nf 1 is as given by eq. (|5[) but now the integration over the apparent energy is from 
11.5 MeV. The time averaged probability for this case is 

(P ee (E u , r, R(t)) = — *— - P ee (E u , r, R{t))dt (10) 
h — h Jti 

where and t\ are determined by the eight bins provided by SK. For each of these bins, 
the SK collaboration has presented the ratio (N^ xp / Nf P98 ) |J. The best-fit values for 
oscillation to sequential neutrinos obtained in this case are 

• Am 2 = 4.26 x lO" 10 eV 2 , sin 2 (2$) = 1.0, X n = 1.12, xLn = 3-96 

The goodness of fit is 55.52%. We find that from the seasonal data almost all of the 
parameter space in the Am 2 — sin 2 (2 , (?) plane is allowed at 90% C.L.. In fig. 4 we 
plot Xmin against one of the three parameters keeping the other two unconstrained 
to illustrate this point. From fig. 4 it is seen that all values of Am 2 in the range 
10~ n — 10~ 9 eV 2 and sin 2 (2 , (9) from very small values to 1.0 are allowed at 90% C.L.. 
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X n below 0.5 and above 3.5 are not allowed at 90% C.L.. Thus, the preliminary seasonal 
variation data does not put any strong constraint on the parameters. For oscillation 
to sterile neutrinos the best-fit values are 

• Am 2 = 4.09 x 10~ 10 eV 2 , sin 2 (2$) = 1.0, X n = 0.99, xLn = 3 - 88 
For 5 (8 - 3) degrees of freedom this solution is allowed at 56.68%. 

7 Conclusions 

In this paper a detailed x 2 -analysis of the 708 day SK solar neutrino data is performed 
assuming vacuum oscillation between two neutrino flavors. This includes 

• A x 2 -analysis using the total rates from the 37 Cl, 7l Ga and Superkamiokande 
experiments. We take into account the experimental and the theoretical errors 
including the correlations among the various theory errors. 

• A x 2 -analysis using the SK spectrum data including the uncorrelated as well as 
the correlated errors among various bins as given by the SK collaboration Q. 
We float the normalisation of the 8 B flux as a free-parameter and determine its 
best-fit value. 

• A global x 2 -analysis of the combined rates and spectrum data treating them 
as independent. For this case also we have allowed the 8 B flux to vary as a 
free-parameter. 

• x 2-aria ly s i s °f the preliminary seasonal data allowing the 8 B flux normalisation 
to vary. Here we consider only the experimental statistical errors. 

We find that the simple two-generation vacuum oscillation scenario can well explain 
the data on total rates and the spectrum but somewhat different best-fit values of Am 2 
are found in the two cases. If one does a global analysis of the rates and the spectrum 
then the goodness of fit is poorer. This is because of the fact that the best-fit value 
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of Am 2 for the spectrum data is one order of magnitude higher as compared to the 
rates value. The sterile neutrino alternative gives a bad fit to the data on total rates 
but for the SK spectrum data as well as for the combined rate and spectrum data it 
gives a marginally better fit than the active neutrino scenario. The sterile neutrino also 
gives a slightly better fit for the preliminary seasonal variation data. Thus we make an 
important observation that if only the SK spectrum or the seasonal data are considered 
then the sterile neutrino actually gives a slightly better fit than the sequential neutrino 
case. It is when the rates in the radiochemical experiments are included that the sterile 
neutrino fit becomes worse. 

Two-generation vacuum oscillation analysis of the 708 days of SK data has also 



been performed in |16|. However their fitting procedure is somewhat different from 
ours. We take into account the total SK rate as well as the spectrum data. We also 
include the theory errors and their correlations in the analysis of the total rates. For 
the SK spectrum data we have taken the bin by bin correlations into account. Since 
they have not included the SK rates we cannot compare our results for the total rates 
with theirs. If we compare our results of only the spectral and seasonal data with 
their analysis, the best-fit values are more or less in agreement though our Xmin are 
somewhat lower. Because the seasonal data does not put strong constraints on the 
parameters, this preliminary data have not been included in our global fit. We have 



not included the data from the measurement of day-night solar neutrino flux |TjJ as 
the vacuum oscillation hypothesis does not generate any day-night asymmetry. 

In order to explain the energy spectrum observed at SK, in the literature the hep 
flux has sometimes been allowed to vary apart from the 8 B flux and it was found 
that a hep flux almost 20-30 times larger than the SSM prediction can well explain 
the high energy part of the spectrum data. In this work we have not considered this 
possibility. Whether one requires such a higher hep flux than the SSM prediction can 



be tested in SK. It was concluded in [[TjJ that even if one allows the hep flux to vary, 
the allowed parameters for the vacuum oscillation hypothesis will not change much. 
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For our global analysis of the rates and the spectrum we have treated these two 
data sets as independent. Since the S B neutrino flux enters the rates as well as the 
spectrum data there can be some possible correlations among these. In a future study 
we plan to examine whether the inclusion of these alters the conclusions by a significant 
amount. 
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Table 1: The ratio of the observed solar neutrino rates to the corresponding BP98 SSM 
predictions used in this analysis. The results are from Refs. || and 0. For Gallium, 
the weighted average of the SAGE and Gallex results has been used. 



Experiment 


Chlorine 


Gallium 


Super Kamiokande 


Observed Rate 


0.33 ±0.029 


0.57 ±0.054 


0.471 ±0.015 


BP98 Prediction 



Figure Captions 

Fig 1. The 90% C.L. allowed region in the Am 2 - sin 2 (2 , (9) plane from the analysis of 
total rates assuming vacuum oscillations to sequential neutrinos. The best-fit point is 
also indicated. 

Fig 2. The 90% C.L. allowed region in the Am 2 - sin 2 (2T?) plane from the SK re- 
coil electron spectrum data for sequential neutrinos (a) and for sterile neutrinos (b). 
The best-fit points are also indicated. 

Fig. 3. The 90% C.L. allowed region in the Am 2 - sin 2 (2$) plane from a global 
analysis of the rates and the spectrum data for sequential neutrinos (a) and for sterile 
neutrinos (b). The best-fit points are also indicated. 

Fig. 4. Variation of Xmin (solid line) (a) with Am 2 keeping sin 2 (2$) and X n un- 
constrained; (b) with sin 2 (2$) keeping Am 2 and X n unconstrained; and (c) with X n 
keeping sin 2 2$ and Am 2 unconstrained. Also shown are the 90% C.L (big-dashed 
line), 95% C.L. (small-dashed line) and 99% C.L. (dotted-line) limits. 
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Fig. 2 
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